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This paper describes a method for computing baseband distortion in 
analog FM communication systems; the method is based on recent theo- 
retical work available in the literature. The input baseband signal is taken 
to be a zero-mean, stationary Gaussian process having an arbitrary power 
spectral density. A variety of graphical results are presented in order to 
demonstrate the utility of this method of computing FM distortion. It is 
shown that the often-used noise loading test does not necessarily represent 
a ivorst-case test. 

I. INTRODUCTION 

Theoretically, analog FM signals generally possess an infinite band- 
width. Thus, when such signals are passed through a linear system 
having a finite bandwidth, some FM distortion must occur. The mea- 
surement of such distortion is costly and very time consuming. Ac- 
cordingly, the development of methods for the computation of FM 
distortion is of practical interest. 

The purpose of this paper is to describe how we used the theoretical 
results derived by A. Mireea, 1 E. Bcdrosian 2 and S. 0. Rice 3 to develop 
a computer program to compute the FM distortion resulting from 
linear time-invariant-filter structures. The input baseband modulation 
is taken to be a zero-mean, stationary Gaussian process having an 
arbitrary baseband power spectral density. 

II. SERIES EXPANSION UNDERLYING THE COMPUTATION 

Consider the analog FM communication system presented in Fig. 
la. The associated mathematical problem for studying FM distortion 
is illustrated in Fig. lb. The problem is to deduce the double-sided 
power spectral density, W e (f), of the output random process, d(t), 
given r(/) and the double-sided power spectral density, W+(f), of the 
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input Gaussian modulation. Once this mathematical problem is solved 
satisfactorily, we can then compute the FM distortion at baseband. 
In the FM case, W*>(f), the power spectral density of <p'(t), is given 
and We>(f), the power spectral density of d'(t), is desired. In the PM 
case, W<,(f) is given and We(f) is desired. However, the two problems 
are closely related because of the following relations : 

WAD =" 2 WM) (i) 

and 

WvU) -«nfi(/). (2) 

In fact, an FM communication system can be designed by using only 
PM equipment, as is illustrated in Fig. la. 

Using Rice's 3 notation, a series expansion of We(f) is given by: 

WeU) = eiAf) + iw*(f)\u(f) + t/*(-/)| 2 

+ - f dpWMWM- p)\T{ P ,f- p) - T*(-p, -/+p)| 2 

+ — f d P f daWMWMWM - p-o) 

24 ./_«, J -„ 

X \S(p, a, f-p-a) +S*(-p, -a, - / + p + <r) | 2 

+ 0(*HF*) (3) 
where 

* = complex conjugate 

dc = dc part of 0(0 

T(p, f- P ) = S( P , f - p) + f daW t (<r)[2S(*, pM~", J-p) 

J —00 

- S(a, f - a) - T(a)T(-<r)S(p, j-p) 

+ S{p + v, /-p-<r)] 

U(f) = r(/) + f* dpW+(p)r(p)S(- P ,f) 

J —HO 

+ f dp! daWMWM{-\Y{p + c)S{-p-<T,f) 

J —00 J —00 

+ r(a)lSS(-a, p)S(-p, f) - S(p, f-p-a) 

+ S(p - ff, f - p)]} 
S(p,a) = T(p + ct) - r(p)r(ir) 
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Fig. 1 — (a) Analog FM communication system. 0'(i) is a distorted version of 
<t>'{t). (b) Associated mathematical model for studying FM distortion. 



and 

S(p, a, v) = T(p + a + v) - V(p + a)T(v) - T(p + v)T(<r) 

- v(a + p)r(p) + 2r(p)r(<r)rOO. 

We shall neglect the dc term, d\ c , since we are mainly interested in 
the continuous part of We(f). In addition, for the rango of parameters 
of interest to us, we have found that the double integral associated 
with the U(f) function can also be neglected. 

Notice that W i> {f l ) = does not imply that W B (fi) = 0. This is 
contrary to the case of a linear system. That is, even if we apply no 
input power in the frequency interval (f lf /i + df), we generally get 
some "intermodulation noise" at the output in this frequency interval. 
Actually, eq. (3) is a truncated form of an infinite series of functionals 
of r(/) and W+(f). However, we shall see that it is possible to select 
system parameters which are of practical importance and which allow 
us to neglect all of the terms represented by O(0 6 TF^,). Accordingly, 
we shall define the signal power S(f)df in the frequency interval 
(f, f + df) at the output to be 



S(f)df = W*U)\U{f) + U*{-f)\Wf. 



(4) 



For the range of parameters which are of practical importance, it 
turns out that S(f) = |W«(/)|r(/) + T*(-/)| 2 . S (/) represents the 
spectral contribution at the output which is free of intermodulation 
noise. We also define the F.M distortion power D(f)df appearing in the 
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output frequency interval (/, / + df) to be 
D(f)df = - I dpWMWM - P ) 

O J —oo 

X \T( P ,f- P ) - T*{- P ,-f + p)\ 2 df 

1 



24J-J P . 



daWMWMW+u- p- °) 

X \S( P ,e,f-p-*)+S*(--p, -«r, -/ + p + <r)|«d/- (5) 

D(f) represents the spectral contribution at the output which results 
from intermodulation noise. 

The two quantities of prime interest in this section arc the signal-to- 
distortion ratio, S(f)/D(f), and the ratio of the average signal power, 
as, to the average distortion power, a 2 D , in the output baseband. The 
latter quantities are defined by 



2 
as = 



<JD = 



S(f)df 



u 2 S(f)df 



D(f)df 



<**D(f)df 



PM Case 
FM Case 
PM Case 
FM Case 



(6) 



(7) 



where B = baseband bandwidth of the modulation. When the ratio 
of a\la 2 D ^ 10, one is usually safe in disregarding the terms labeled 
O(0 8 W ) in eq. (3).* 

III. NUMERICAL METHODS EMPLOYED 

An input power spectral density, WV(/), which is often used when 
measuring FM distortion is the uniform spectrum, given by 



(2ttZ)) : 



WAf) = { 2B 



' 0, 



l/l ^B 
f\>B 



(8) 



f Equation (3) is a special case of a much more general equation which was recently 
reported by E. Bedrosian andS. 0. Rice in the Proc. IEEE, 59, No. 12, pp. 1688-1707, 
eq. (14) and Section IVc, December 1971. 
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where D = RMS frequency deviation and B = baseband bandwidth. 
From eq. (1), we have 



D- 



S\^B 



WAf) =\2BF' ~ (9) 

I- , \f\> B 

When such a uniform W+'(f) is used to measure FM distortion, the 
measurement is referred to as a "noise loading test." The noise loading 
test is used, for example, to estimate the FM distortion in microwave 
relay systems resulting when thousands of telephone channels are 
multiplexed to form a composite baseband signal. 

A bandlimited form of Wf(f) is very convenient numerically, since 
it serves to convert the infinite limits of integration in eqs. (4) and (5) 
into finite limits of integration. However, if we attempt to evaluate 
equations (4) and (5) using a bandlimited W*(f) such as is given in 
eq. (9), we would run into difficulty whenever the argument of W^(-) 
is equal to zero. In order to circumvent this apparent difficulty, we 
have selected an integration grid such that the argument of W^,{-) is 
never allowed to be zero. Equations (4) and (5) are then numerically 
evaluated by using a combination of Simpson's rule and the trapezoidal 
rule. 

The particular integration grid used was obtained by setting 

p = (2i + 1)A (10) 

a = (21 + 1)A (11) 

/ = (2tj + 1)A (12) 

A = (2k + 1)-' (13) 

B = 1.0 (14) 

where i, I, n, k are integers. In most of our numerical evaluations, 
k = 20. 

IV. NUMERICAL RESULTS 

4.1 Test Cases 

In order to test the operation of the computer program, we evalu- 
ated D(f) for the case when TIV(/) is uniform and r(/) is the transfer 
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Fig. 2— Baseband intermodulation spectrum resulting from FM distortion. 



function of a simple RC filter. That is, W*(f) is defined by eq. (9) and 
r(/) is denned as 

r(/) = X — (is) 

where 

fc = 1.0 

D = 1.0 

B = 1.0. 

Figure 2 shows the resulting computer plot of oj 2 D(/). Figure 3 shows 
a plot of 10 log [S (/)/£(/)] for / = 0.084B, 0.36B, and B, as a func- 
tion of the RMS frequency deviation D. These results compare well 
with both the theoretical and experimental results which are presented 
in Refs. 2, 4, and 5. 

As a final test case, we present the results for the case when W* (/) 
is still defined by eq. (9), but T(f) now represents a 3-pole Butterworth 
filter with some mistuning, in which case 



r(/) = 



1 + 2(tf ) + 2(tf,)» + (tg )' 
1 + 2(tf) + 2{i$Y + (t*)' 



(16) 
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Fig. 3 — Signal-to-distortion ratio at particular frequencies resulting from T(/) 
and HV(/)- 
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Fig. 4 — Signal-to-distortion ratio resulting from T(/) and U>.(/) for particular 
values of frequency offset or mistiming. The points are from the theoretical approxi- 
mation given as eqs. (17) and (18). 
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Fig. 5 — Signal-to-distortion ratio resulting from T(/) and Wp{$). 



a = 



/-/. 



/. 



£o= - 



/. 



1.0 



B = 1.0 

D = 3.12 

/. = 17.0 

/. = 0, 1, 2. 

Figure 4 presents the computer plot for this case. The results compare 
very well with experimental and Monte Carlo results presented in 
Ref. 4. 

A theoretical approximation for the above case, with ^ | /| ^ B, 
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Fig. 6 — Signal-to-distortion ratio resulting from r(/) and two forms of HV(/). 
was developed by Rice 3 and is given by 



0£>- 



D l 



and 



+ — (3£ 2 -/ 2 )(a 3 ,) 2 (17) 



8(f) = W,(f) (18) 

where X ni is the imaginary part of X„ and \ n /n\ is the coefficient of 
/" in the power series expansion of In r (/) 

x 2l = 2/ e /; 8 

X 3 , = -2/; 3 

\u = 48/ e /7 5 . 
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Some points which were computed from this theoretical approximation 
are indicated in Fig. 4. 

Having verified the soundness of the computer program with the 
foregoing test cases, let us present some new results. 

4.2 New Results 

In this section, we shall present some new results which were ob- 
tained by using the above methods. These results will also help to 
demonstrate the general types of FM distortion problems that can be 
analyzed. 

4.2.1 n-pole Chebyshev filter 
Let W*>(f) be uniform as given by eq. (8) with T(f) given by 

n 

n too - 8 k ) 

r(/) = ^ (19) 

fl {H - 8k) 
where 
St =- S1 n[(2fc-l)|Jsinh[isinh-'Q] 

+ icos|(2fc - 1)— cosh -sinhrU-J , * - 1, ■ • • , n 

B = 1 = baseband bandwidth 
2/ c = K2B(l + m) = K times Carson's rule 
= filter bandwidth 

4D . . , 

m = — = FM modulation index 
B 

R = 10 log (1 + b 2 ) = in-band ripple 

f e = offset frequency or mistuning. 

r(/), defined by eq. (19), represents an n-pole Chebyshev filter. Some 
results for this case are presented in Figs. 5, 6 and 7. Also, a computer 



COMPUTING BASEBAND DISTORTION 



637 





14 














\ 1(0 = 


9-POLE 


CHEBYSHEV FILTER 




40 


- \ f c = 


(2/3) (16/15) 


B(1 +m) 






\ ° = 


RMS FREQUENCY DEVIATION 






\ \ B = 


BASEBAND 


BANDWIDTH = 1.0 




36 


~ \ \ m = 


4D/B = 


FM 


MODULATION INDEX = 2.0 






\ \ R = 


0.3-dB 


RIPPLE 




32 


\ \ f e = 


n 






Z 












O 

< 

a. 


28 


~~ \ \ 








7 

o 


24 










h- 
o 












5 


?0 








^^W„y(f) =UNIFORM SPECTRUM 


i 

o 
- 

i 

_j 
< 
z 


16 
12 








^"^- k _____ — 


— — =21.0 dB 




B 








. W^'(f) = RC SPECTRUM S^ .2 
"\f = 0.1B^/ -i--1B.7dB 
D 




4 


1 1 1 






I I I I I I 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

NORMALIZED FREQUENCY, f / B 

Fig. 7 — Signal-to-distortion ratio resulting from r(f) and two forms of HV(/). 



plot is presented in Figs. G and 7 for the case when W r (f) is an RC 
spectrum 1 given by 



WAS) = « 



CZtrD)- 



B 

2/ tan- 1 { — 
/o 



-3T 



where 



/o = 3 dB bandwidth. 



/I -£ 5 



/I >B 



(20) 



An RC spectrum is often used to model a video signal. Notice that, as 
fo — * °° , the RC spectrum approaches the uniform spectrum as given 



r II '*'(/) is the spectrum produced by passing bandlimited "white" noise through 
an RC filter. 
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Fig. 8 — Signal-to-distortion ratio resulting from r(/) and 11 V(/). 



by eq. (S). From Figs. 6 and 7, we see that more FM distortion results 
when W*> (/) is an RC or video spectrum than when W> (/) is uniform. 

4.2.2. n-pole, equalized Chebyshev filter 

If the phase function associated with eq. (19) is taken to be zero/ 
T(/) can be written as 



r(/) = 



- 1 + WTUSoW 

L i + bm(o J 

where T n (£) is a Chebyshev polynomial given by 

[cos [n cos -1 (£)], |f | ^ 1 



(21) 



TnU) = 



cosh [n cosh l (£)], 



*l > 1 



T(/), defined by eq. (21), represents an ??-pole, equalized Chebyshev 
filter. Some results for this case are presented in Fig. 8. By comparing 
Figs. 5 and 8, we can determine the effect of equalization on FM distor- 
tion. In this case, equalization does not reduce the FM distortion 
significantly. 



* Or linear in frequency since time delay is unimportant here. 



COMPUTING BASEBAND DISTORTION 



639 




0.3 0.4 0.5 0.6 0.7 

NORMALIZED FREQUENCY, f /B 



Fig. 9 — Signal-to-distortion ratio resulting from r(/) and TJV(/). 

4.2.3 7i-pole Butterworth filter 
Let Wf(f) be uniform as given by eq. (8) with T(f) given by 



r(/) = 



n (po 



(22) 



where 



II (w + Vk) 

fc-1 

V k = (2r/J exp ft j ( — ^ 1 X\ , k - 1, 2, 



5=1 = baseband bandwidth 
2/« - A'25(l + m) 
4D 



m = 



B 
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Fig. 10 — Signal-to-distortion ratio resulting from r(/) and UV(/) . 

T(/), defined by eq. (22), represents an w-pole, Butterworth filter. 
Some results for this case are presented in Fig. 9. By comparing Figs. 
5 and 9, we see that the 3-pole Butterworth filter produces much less 
FjVI distortion than does the 9-pole Chebyshev filter. 

4.2.4 n-pole, equalized Butterworth filter 

If the phase function associated with eq. (22) is taken to be zero, 
or linear in frequency, T(f) can be written as 

r(/) = 



l + 



m 



(23) 



r(/), defined by eq. (23), represents an w-pole, equalized Butterworth 
filter. Some results for this case are presented in Fig. 10. By comparing 
Figs. 9 and 10, we see that equalization reduces the FM distortion 
significantly in this case. 
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4.2.5 Echo, envelope delay, Butter worth filter 
Let W+> (/) be uniform as given by eq. (8) with T (/) given by 
1 + re'" 7 "! 
+ r _ 
(echo) (envelope delay) 

1 +2(tf ) r-2(t{o) 8 + (ih) 9 ' 



nf) - 



K 



[exp W& 2 co 2 + 6 3 co 3 )}] 



X 



where 



_ 1 + 2(t£) + 2(tfl* + (t{) 
(3-pole Butterworth filter) 

/• = amplitude of echo 
r = time delay of echo 
bi = linear envelope delay constant 
6 3 = quadratic envelope delay constant 



it) 3 J 



* = 



fc 



(24) 



/ f = frequency offset or mistuning 
D = RMS frequency deviation 
B = 1 = baseband bandwidth 
2/ c = if2fi(l + m) — X times Carson's rule 
= filter bandwidth 

4D 

m = — = FM modulation index. 
B 

Some results for this case are presented in Fig. 1 1 . 

Results are presented in Fig. 12 for the case when r(/) is given by 
eq. (24) and W+(f), rather than W *•(/), is uniform and given by 

r3D 2 

W*(f) = JiS*' l/l ~ B (25) 

[ , l/l > A 

This is the case of a noise loading test applied to a phase modulated 
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Fig. 11 — Signal-to-diatortion ratio resulting from r(/) and WyU). 



system. From eq. (1) we have 

|-(2tZ)) 2 3/ 2 
WAf) = 






2B 3 





f\SB 

f\ > B. 



(26) 



In this case, W*> (/) peaks at / = B in contrast to the RC spectrum 
given by eq. (20), which peaks at / = 0. These results lead us to an 
interesting question. Given W r (f) and r(/), can we choose a predis- 
tortion characteristic such that the shape of S(f)/D(f) is most suit- 
able for a particular communication system? However, we have not 
investigated this question. 

We can compare the distortion resulting from PM and FM systems 
by comparing Figs. 11 and 12. In fact, if the results in Fig. 12 apply to 
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Fig. 12 — Signal-to-distortion ratio resulting from r(/) and W+(f). 

a Phase modulation system, then the results in Fig. 11 apply to the 
corresponding Frequency modulation system. 



4.2.6 Filter characteristic determined experimentally 

In all of the results presented above, the filter characteristic T(f) 
was specified mathematically. However, many situations arise for 
which the measured amplitude and envelope delay characteristics are 
available in graphic form. Of course, one may try to fit a suitable ana- 
lytical expression to these experimental points and proceed as above. 
However, there is no need to develop such an analytical expression. 
As our final case, we shall present some results which illustrate this 
point. 

Let the amplitude and phase of T (/) be given by the experimental 
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Fig. 13 — Experimental amplitude and phase characteristic specifying r(/). 
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Fig. 14 — Signal-to-distortion ratio resulting from V(f) and W r (f)- 
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rif) = AMPLITUDE AND PHASE OF FIGURE 13 
W,//(f) = RC SPECTRUM, f o = 0.1B 

D = RMS FREQUENCY DEVIATION 
B=BASEBAND BANDWIDTH = 1.0 
m = 4D/B » 1.0 
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Fig. 15 — Signal-to-distortion ratio resulting from r(/) and WV(/). 

curves shown in Fig. 13. Lot Wf(J) again be uniform and be given by 
eq. (8). Figure 14 presents some results for this case. 

Now let W^if) be an RC spectrum given by eq. (20). Figure 15 
presents some results for this case. We see that considerably more dis- 
tortion is indicated in Fig. 15 as compared with Fig. 14. 

Thus, the often-used noise loading test which makes use of a uniform 
spectrum may not represent a worst-case situation as far as FM dis- 
tortion is concerned. A theoretical proof of this interesting fact is pre- 
sented in the appendix. 

V. CONCLUSIONS 

Equations (4) through (7) , together with a digital computer, can be 
used to compute the FM distortion resulting from passing FM waves 
through linear networks. To demonstrate the utility of the program, 
we have presented a variety of results in graphic form. 

From this work, it is apparent that the often-used noise loading test 
does not necessarily represent a worst-case test. This was demonstrated 
for a system in which the modulating signal is a video signal. It is also 
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apparent that predistortion may be useful in reducing FM distortion. 
This is in contrast to the use of conventional pre-emphasis, which is 
applied to combat RF noise. 

Some important sources of FM distortion which were neglected in 
our analysis are AM-to-PM conversion, baseband and RF noise, and 
adjacent channel interference. 

VI. ACKNOWLEDGMENTS 

It gives me great pleasure to acknowledge helpful discussions with 
S. 0. Rice. I also wish to express my thanks to J. N. Lahti and 
A. Yagoda for some useful comments, and to Miss R. Wright and F. L. 
Crawford for programming the digital computer. 

APPENDIX 

Theoretical Comparison of FM Distortion Resulting from Video and 
Uniform Spectra 

The purpose of this appendix is to present a theoretical argument 
which shows that a video spectrum can lead to more FM distortion 
than a uniform spectrum. That is, a noise loading test which makes use 
of a uniform spectrum may not represent a worst-case test when the 
information source is a video signal. 

Let the video signal <f>'{t) be represented as a zero-mean, stationary 
Gaussian process having power spectral density W^if) given by 



WAfl - i 



2/o tan 



■0 







f\SB 






(27) 


f\>B 





where 



D = RMS frequency deviation 
B = baseband bandwidth 
/o = 3-dB bandwidth. 

Let r (/) be given by 

r(/) = exp [i(6 2 o) 2 + 6 4 w 4 )] (28) 

where 

b 2 = small linear envelope delay constant 

b t = small cubic envelope delay constant. 
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By applying eq. (24) of Rice, 3 the leading term of D(/) can be ex- 
pressed as 

D(f) = 2->(X 2l + 2-i/)%,) 2 ( dpWMW+U - p)p*if ~ p) 2 (29) 

J —oo 

where 

X«- (2!)(2,r)'6 s 

X 4 , = (4!)(2x) 4 6 4 . 

By taking S(f) = W*if) = W r (f)/e> 2 and evaluating the integral for 
D(f), we find that 



RU) = 



x 



F 



where 
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(30) 



R(f) represents the distortion-to-signal ratio for a video spectrum 
divided by the distortion-to-signal ratio for a uniform spectrum. If we 
can show that R if) > 1 for particular values of /o, the 3-dB bandwidth 
of Wfif), then we can conclude that a video spectrum can produce 
more FM distortion than a uniform spectrum. 

In order to show that R(f) can be greater than unity, consider the 
important frequency rango < \F\ < 1. For this baseband frequency 
range, eq. (30) yields 

4 

lim Rif) = — > 1. (31) 

F -o 2 — |F| 

A plot of Rif) for various values of F is shown in Fig. 1G. 
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r(f)-EXP[i (b 2 eu 2 + b 4 a) 4 )] 
W0'(f)= VIDEO SPECTRUM, f = 3-dB BANDWIDTH 
B = BASEBAND BANDWIDTH 
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NORMALIZED FREQUENCY, f/B 
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Fig. 16 — R (/) denotes the distortion-to-signal ratio for a video spectrum divided 
by distortion-to-signal ratio for a uniform spectrum. 

Accordingly, we conclude that a video spectrum can produce more 
FM distortion in the baseband frequency range than a uniform 
spectrum. 
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